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Abstract. Dynamical invariants and statistical properties of the quantized 
electromagnetic field in nonstationary linear media (dielectric and/or conductive) are 
£^ ' considered in the framework of Choi-Yeon quantization scheme. It is shown that in 

the eigenstates of linear dynamical invariant the Robertson-Schroedinger uncertainty 
^h' relation is minimized, both for the photon annihilation operator quadratures and for 

the electric and magnetic field components. The time evolution of initial Glauber 
coherent states and Fock (photon number) states is considered. On an initial coherent 
■ state the medium conductivity and the time-dependent electric permeability both are 

CO , shown to act as squeezing and correlating factors. 

in 
in 

On 

O 1 1- Introduction 

O ' 

Although the remarkable achievements have been made in classical and quantum 
optics, yet probably many properties of light still remain to be uncovered. For deeper 
understanding of its nature and more precise manipulation the light needs be quantized. 
The method of quantizing light propagating in free space is well known and appears in 
most quantum optics text books (see e.g. the book of Louisell [lj, Scully and Zubairy 
[2], Walls and Milburn [3j). Practically the same is the light quantization procedure 
in a stationary homogeneous and isotropic dielectric media, the only new feature being 
the reduced light velocity. The quantization of damped light is somewhat more subtle. 
It involves the quantum description of a single (or several) field oscillator (s) interacting 
with a reservoir with a large (infinite) number of degrees of freedom (see e.g. pi] and 
references therein) that makes the calculations rather lengthy. Besides, the Hamiltonians 
they used in the development of the theory are somewhat assumed ones, rather than 
having been derived in a consistent way from the classical electrodynamics. Recently 
publications by Choi and Yeon appeared [5j [6J |7J EJ [9] (see also |10[ [TT]) where a scheme 
was proposed for quantizing the damped light in conducting (and nonstationary) linear 
media without link to a reservoir, resorting to the Caldirola-Kanai Hamiltonian [12j and 
to the Lewis-Riesenfeld dynamical invariant theory 
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Our main purpose here is to construct, in the framework of the Choi-Yeon 
quantization scheme [TJ [8], Robertson-Schrodinger correlated intelligent states [TH [15] 
of electromagnetic field (EMF) in nonstationary linear media and to examine some of 
their properties. We also consider the time evolution in such media of initial Glauber 
coherent states (CS) and Fock states and show that the conductivity and/or the electric 
permeability act as a correlating and squeezing factors. We resort to the linear dynamical 
invariant method, developed in papers by Malkin, Man'ko and Trifonov |16[ ITT] (see also 



2. Quantization of EMF in linear media 

As in papers [TJ El [9] we consider EMF in nonstationary (isotropic and dispersionless) 
linear media with no charge sources (the case of stationary conductive linear media was 
considered earlier in [5j[6]), 

B(r, t) = n(t)H(r, t), D(r, t) = e(t)E(r, t), j = a(t)E, (1) 

divD = 0, divB = 0, (2) 

d 

rotH = —D + a(t)E, (3) 

votE = —ttB. (4) 

dt v 1 

We fix the Coulomb gauge with vanishing scalar potential, 

B = rotA, E = (5) 

dt v 1 

and write the equation for A(r, t) ( the term e being omitted in papers [7J [S]) 

V 2 A-v(a + e)—-e^—=0. (6) 

As usual (see e.g. the books [TJ El |3]) we expand vector potential A(r,t) in terms of 
mode functions Ui(r) = ei^ui^(r) 

A(r,t) =J2 e ht u i,d r M t ), ( 7 ) 

where, for further convenience, we introduced the polarization vectors e^. The mode 
functions W/,^(r) are supposed to be eigenf unctions of the Laplace operator [IJ[2l|3] (ojoj 
are constants of dimension of frequency), 

V 2 + ^) u«(r) = 0, (8) 

to form a complete orthonormal set and Ui to satisfy the transversality condition, 
V ■ U\ =0. For example, in case of EMF in a cavity of the form of rectangular cube of side L = y 1 / 3 
the polarization mode in xi-direction the appropriate mode is Ui(r) — e x sm(l y i:y/L) sin(l z Trz/ L), 
where l y , l z are integers. In the case of periodic boundary conditions, the mode functions may be taken 

as mci! 

un(r) = V- 1/2 e^ exp(±ik l ■ r), (9) 
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where ki is the wave vector, fc/ • e;^ = 0. In three dimensions ki — 2ir'^2 i h e i, i — 1, 2, 3 where li are 
integers and ej are orthogonal unit vectors. The set 1%, la, gives a mode of a given polarization. 

From Maxwell equations (j2J) - (jlj) it follows that (in case of linear media ([I])) the 
time-dependent factors qi are to obey the following linear equation (furthermore, unless 
otherwise stated, we suppress the polarization index £) 

d 2 qi WjA % 



where (the term e being omitted in [3 [8] ) 



9 ' 

ere 



r(t) 



a + e 



(10) 



The frequencies and the damping parameter V are time dependent if the medium 
parameters are. Noting that the eq. ffTOl) could be obtained from the Hamilton function 



H, 



fc,f 



A(t) 



£o 

r(t')d*', 



the authors [TJ [8] perform the quantization introducing the Hamiltonian operator 



H 



E 



2e yi 



eoufit) 



(12) 
(13) 

(14) 



were p;/ are Hermitian operators satisfying canonical commutation relations, 

[qt,Pi>] = ihSw, [q h qi'\ = \PhPl') = 0. (15) 

For constant A and u>i the one mode Hamiltonian (and the Hamilton function (TT2l ) 
can be recognized as the reknown Caldirola-Kanai Hamiltonian [12] . considered for 
description of damped mass oscillator motion. The vector potential operator reads 

A(r,t) = J2 u i(r)qi. (16) 



Replacing it in fl5]) one obtains the quantized electric and magnetic fields 
1 



E(r,t) 
B(r,t) - < 



e- A ^J2Mr)Pi 



1 



x urfrjqi. 



(17) 
(18) 



Using the time derivatives of operators qi, pi of the form 

ddi i ~ 

!t = -h [ohHl > 

we check that all Maxwell equations (J2])- 
eE, H and B = /iH. 



(19) 



are satisfied by operator fields E, D 
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It worths noting here that the Hamiltonian (1141) (the Hamilton function (TT2")) ) that 
governs the time-evolution does not coincide to the field energy operator W (the field 
energy W), the latter being given by the expression 
i / -2A(t) \ 

w=^< t n—^pi+^n^J2^- ( 2o ) 

1 i V £ o / I 

This difference is a feature of dynamics of open systems. For the case of classical damped 
harmonic oscillator with constant T and u it was first noted by Bateman in 1931 [19J. 
The stationary damped oscillator was later quantized by Caldirola and Kanai |T2] using 
Bateman classical Lagrangian [19j. Now we see that quantized EMF in conducting 
media could be regarded as a set of (noninteracting) damped oscillators. 



3. Dynamical invariants and Ermakov equations 

Dynamical invariants I of quantum system with Hamiltonian H are defined as solutions 
to the equation 

§4[/,ff]=0. (21) 




. If, . a 

ep H — eo — e e 

a \ 2a 



The canonical commutation relations (1151) show that quadratic in q and p Hamiltonians 
admit linear in q and p dynamical invariants. In refs. [T7] a family of (non-Hermitian) 
invariants A for the general nonstationary quadratic Hamiltonian 

H = ~ [a(t)p 2 + b(t)(pq + qp) + c(t)q 2 ] , (22) 

have been constructed in the form 

/ " ^ 1 (23) 

where e is any solution of the second order equation (classical oscillator equation) 

e + fi 2 (t)e = 0, (24) 

n 2 = ac + b ^ + ^__^_ b 2_ k (25) 
a 2a 4cr 

The commutator [A, A^\ reads 

[A,ti]= l -{ee*-e*e) = l -w. (26) 

This shows that if the Wronskian w of the equation (|24|) is fixed to —2i (as in [TBI ITT]), 
then the invariants A and are boson ladder operators. (It is worth noting that real 
solutions e(t) are also admissible, they correspond to Hermitian invariants). The relation 
w = —2i is identically satisfied with e of the form e = |e| exp(i J dt'/\e(t') | 2 ). Now eq. 
(|24|) leads to the Ermakov equation [20J for \e(t)\, 



d 2 \e\ 
dt 2 



+ n 2 (t)\e\---, = 0. (27) 
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At 



a(t) = e^e-W, b(t) = 0, c{t) = e co 2 (t)e A ^ (28) 

the Hamiltonian (122]) recovers the nonstationary one-mode Caldirola-Kanai Hamiltonian 
Hi, eq. pijl . Therefore at ( 128]) the ladder operator linear invariants A(t), eq. ( 123]) are 
dynamical invariants for the damped oscillator. Thus for each mode / of the quantized 
EMF we have dynamical invariants in the form of boson annihilation operators 



.4, 



-±A(t) 



em - e e A ^ 



(29) 



^/2hF 

where e; satisfy the equations 

€i + Qf(t)ei = 0, (30) 

n?(t) = - 1a - ^A 2 , (31) 

and £/e* — e*e; = — 2z. Note that under the latter conditions all ||q| are solutions to the 
Ermakov equaton (l27j) . 

For general quadratic Hamiltonian ( 1221) the linear invariants A(t) and their 
Hermitian combination A\t)A(t) have been diagonalized in |17| : 

A(t) \a; t) = a\a; t), A\t)A{t)\n-t) = n\n;t), (32) 

where a G C and n = 0,l,2,.... For our quadratic Hamiltonians Hi and the invariants 
Ai, A\Ai formulas (44) and (46) of the first paper in |17| produce the eigenf unctions 

ip ai (qi,t) = (qi\ai;t), ip ni {qi,t) = (qi\rif,t) (a = e 1 e~ A ® = a(t)), 

~2 a t 



i> ai (qi,t) = i/j (qi,t)exp 




ah ei 



qi 



— a, 

2e t 1 2 



\ni/2 



"ni< 



H n ,(xi) 



xi 



qi 



(33) 



(34) 



where H n (x) are Hermite polynomials, ^o(qi,t) are the ground state wave functions, 



ipo(qi,t) = [ei(nah)* 



cxp 



2ah \ei 2a 



Qi 



(35) 



These time- dependent wave functions are normalized solutions to the Schrodinger 
equation ihd t ip = Hi%p with Hamiltonian Hi given by eq. (j!4p . Since A(t) and A\(t)Ai(t) 
are dynamical invariant, the eigenvalues ai and n\ are constant in time. 

The system of \cxi',t) is overcomplete in the one mode Hilbert space Hi (the set of 
\ni;t) being complete): 



- [ \ai;t)(t;ai\d 2 ai = ^ \m;t)(t;ni\ = 1 
71 J ni 



(36) 



According to the terminology of Refs. |16[ [T7] the states \af,t) may be called CS of 
nonstationary system with Hamiltonian Hi, eq. ( 1T4]) . Since the EMF Hamiltonian H, 
eq. (TH]) . is a sum over /, CS for EMF with finite number of modes are product over I 
of one mode CS \af,t). 
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Remark 1. To make connection with the Choi and Yeon treatment [71 [S] we put 



e,(t) = 



1 



A 



2M t 
M, 



=>»TJ(*) 



(37) 
(38) 



Then we find that the two real parametric functions 7; (t) and M; (t) have to obey the following system 
of (second order nonlinear) equations 



M, 




(<tf-w?(t))-Afi 



(T + £ 



(39) 



(40) 



Up to the term Mii/e these equations do coincide with eqs. (17), (18) in Ref. [7] and eqs. (7), (8) in [5]. 
Evidently this term was omitted in [7J [5] • Under the substitutions (f3"Tf , (|38[) our linear and quadratic 



invariants Ai(t), A\(t)Ai{t) and their eigenstates \ai;t), \ni\t) will coincide with Choi- Yeon operators 
exp(i , Yi(t))ai(t), a»(i)aj(t) and their eigenstates \ai(t)), \<fi n ,i) H[E]> if the parametric functions Ji(t), 
Mi(t) in [7, 8j are subject to the equations (|39|) . (j40|) and their T(t) is replaced with our T(t), eq. (jTTJ) 



4. Evolution of CS and photon number states in linear media 



Let a ojQl be photon annihilation operator corresponding to the frequency ojqj- In terms 
of Hermitian operators qi and pi its expression reads 



^0^0,/ 

2h 



qi + 



-Pi 



(41) 



One has [a Wo , 



a 



,1' ^0,iJ 



1. Glauber coherent states (CS) |cfy) of one-mode EMF are defined 
as eigenstates of the photon annihilation operator a WQ , 

o^K) = ->■ K) = ex P ("j^wo,; _ a ^i ,i) I )' ( 42 ) 

CS of the iV-mode field are given as product of the /-one-mode CS. 

The ni photon states \n{) (the Fock states or the number states) are eigenstates of 
photon number operators hi = oi^ a Wot , 



ni\ni) = ni\ni) -» \ni) 



'nil 

The normalized wave functions of these states are 

i'aMl) = ^o(*) eX P 



-|0). 



<*iqi 



4>n,(qi) = Mil) 



'nil 



Hniiqiy/UojEo), 



OF 



(43) 



(44) 



(45) 



where H n (x) are Hermite polynomials, and i[>o(qi) is the ground state wave function 

1 



exp 



^O^CM 2 

9/ 



2h 



(46) 
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The photon statistics in Glauber CS \ai) is Poissonian with mean (aj|nj|) = |cfy| 2 . 
They are remarkable also with minimal fluctuations in q\ and pf. the variances of the 
dimensionless quadratures of photon ladder operators Qi = yJ eou j° ,l qi, Pi = \J £Q ^ t n Pi 
are equal at the lowest possible level, 

(AQ,)1 = m)l t = ~, (47) 
minimizing the Heisenberg uncertainty relation (UR) 



1 



4 



2 1 



(AQi) (APi) > — ([Q h Pi]) =-. (48) 



4 



In CS \oti) the %-pi covariance Cov(qi,pi) ai is vanishing. 

Glauber CS and number states are both time-stable in vacuum and in stationary 
nonconductive linear media as well. In these cases the fluctuations of electric and 
magnetic field in CS are also minimal, as one can see from eqs. ( I17jl . (iTBl . In conductive 
and/or nonstationary linear media Glauber CS and photon number states can be shown 
to be both time-unstable. Our next aim is to examine their time-evolution. 

Our time-dependent states \oti;t) and \ni;t) are completely determined by a solution 
ei(t) to auxiliary classical equation (130]) . Therefore initial conditions q(0), e(0) of e(t) 
determines uniquely the evolution of quantum states of EMF in nonstationary media 
(Tjrj. One can check that at 

q(0) = — I— e l (0)=t^W) (49) 
V "/(U) 



the wave functions i/j ai (qi,0) and "0 n; (g/,O) recover the Glauber CS ip ai (qi), eq. 
and Fock states if} nt (qi), eq. ( j45l) with uqj = 17/(0). (If in addition A(0) = then 
Qi(0) = u>i(0).) This means that the initial CS \ai) and the initial photon number states 
\ni) in linear media ([1]) evolve into states \ai;t) and \nf,t). At t > the states \ai;t) 
and \nf,t) deviate from the form of Glauber CS and photon number states, that is they 
are no more eigenstates of the Schrodinger operators a W0l and a u)Ql correspondingly. 
They remain eigenstates of the dynamical invariants Ai(t) and Aj(t)Ai(t) instead. As a 
result the photon statistics in \af, t) is no more Poissonian, the fluctuations of Qi and Pi 
deviate from their minimal value of 1/2 and do not minimize anymore the Heisenberg 
UR. It is worth noting at this point that the invariant A/(i) is a non-Hermitian linear 
combination of the photon creation and annihilation operators a LOol and o£, o; - 

The properties of states with Gaussian wave functions like ip ai (qi, t) (eigenfunctions 
of non-Hermitian linear combinations of photon creation and annihilation operators) 
are well examined in the literature under the names: CS of nonstationary systems 
[THldT], Stoler states [2T] . squeezed states (SS) [22], two-photon states [23J, correlated 
states [14] . generalized intelligent states [15J. The name SS refers to their property 
of "squeezing" the fluctuations of Q or P. The name 'correlated CS' underlines the 
existence of nonvanishing q-p covariance and the name 'generalized intelligent' stresses 
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on their 'generalized minimum uncertainty property' in the sense of minimization of the 
more general UR of Robertson-Schrodinger 

h 2 



(A^) 2 (A^ 



(50) 



In the time-evolved CS \ai\t) we have [25 

ha 

Pi 



"7 



2 

h 

2a 



1 

p! 



Pi(t) 



2a 



Pi 



(51) 



and the covariance [26] (let us remind that now a = e 



e / 



Cov(ft,pj) a| 



-2 Pl \ Pl + Ya Pl 



(52) 



It is seen that these three second moments do minimize UR ( joOjl . In the Yuen [23] 
(m, ^-parameters Ai = ui{t)a Uo , +'Sji(t)a^ , the variances of the dimensionless canonical 



operators take the form 



(AQ* 



«7 



1 



2' w Wl v /ai 2' 
which clearly show that the fluctuations in Qi (Pi) can be 'squeezed' much below their 
ground state value of 1/2 when vi —> U\ [v\ — » —U\). 

In the case of stationary media ( constant e, /i and a) the 'frequencies' Qi are 



constant, Qf 



LOT 



a 2 /4e 2 . In this case the solutions to eq. (130]) are 



ej(t) = exp(zn;t), pi = Q l = constant, (53) 

and our formulas (15 lj) . (l52l) show that the g;-p; covariances are time- independent, 
Cov(qi,pi) a = ha/AeQi, and 



(AgOi(Ap 



/a; 




(53a) 



If in addition a = (i.e. A = 0), then the g;-]3« covariances are vanishing and the 
fluctuations of the dimensionless Qi and Pi become equal to 1/ \[2. 

Thus on the initial Glauber CS of EMF the medium conductivity and the time 
dependence of electric and/or magnetic permeability may act as qrpi correlating and 
squeezing factors. For such stationary conductive media note the dumping factor 
exp(— at/e) in the variances of qi and the amplifying factor exp(o~t/e) in the variances 
of pY 

Remark 2. Exact solutions of eq. (|24l) for e(t) are known for several cases of f2(t). For example 



if n(t) decreases in time as f2(£) = 



t) the general solution is 



(54) 

)\ cos(2i) the eq. 



e(t) — ciy/r(t) cos(s lnr(£)) + c 2 \Jr(t) sin(s lnr(£)), 

where r(£) = £+l/o;o, s = \/3/2, and c 12 are arbitrary constants. For f2 2 = a. 
is known as Mathiew equation. Let us note that the new auxiliary functions e;(£) := ej(£) exp(— A/2) 
have to obey the classical damped oscillator equation (TlUl) : 

l[ + ke[ +u}f{t)e' = 0. (55) 



DIRSCS in NLM- 



9 



For A = r], f\ = and constant u>i solution for ej(i) in terms of Bessel functions J 1; Y\ is also 
known [10J . Noting the analogy of eq. (|24p to the one dimensional stationary Schrodinger equation 
(p/2m + U(x))4>(x) = Eip(x) one can find solutions e(t) in all cases of U(x) where is known, using 
the correspondence fl 2 (x) = 2m(E — U(x))/H |16) . 



5. Statistical properties of quantized field 



For the study of statistical properties of quantized field in the time-evolved CS \cti]t) 
and number states \rif,t) it is convenient to express Hermitian operators qi, p\ in terms 
of the dynamical invariants Ai (t) , A\ (t) . We have 



(56) 



where 



vi = m 



-A(*)/2 



2e 



1/2 



(I- 



-H¥) 1/2 U-^ 



(57) 



We shall consider the case of periodic boundary conditions with complex mode 
functions (Q. With these modes the Hermitian operator of the vector potential, which 
obeys the equation Q) takes the form 



A(r,t) 



2e 



exp 



u. 



^(r)eiAi(t) + h.c. 



(58) 



where is the polarization vector of mode I. The operators E and B are obtained 
from this A via the relations ([5]): 



E(r,t) = y/^e-l*® E u e Ii€ (|Ac, - e,) «&(r)4(f) + 
B(r,t) 



ft_ p -±A(t) 
2e 



/i/^-' x r '/.5 w^(r')e;A z (t) - /i.e. 



(59) 



The commutators between the j and m components of Ei(r,t) and Bi(,t) are C- 
numbers, vanishing for j = m: 



[£^r,i),A, m (r,t)] 



. h 



(60) 



J2 e i,u( k i x e ^)i Re(e /e ; - |e,| 2 A/2)<5 im . 

One mode field. We shall calculate the first and second statistical moments 
of electric and magnetic field in the time-evolved CS \otut) in the case of one mode 
polarized field propagating along x-direction. In terms of parametric functions 7z(t), 
Mi(t) (see Remark 1.) the first moments of E\ and Bi and the means (Ef) a , {Bf) a have 
been calculated by Choi [8J. 
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Using §§§ and the eigenvalue property 

Ai(t)\ai;t) = ai\ar,t), 
we get the first moments in the forms 

(El) ai = 



10 



h 



2e a V 



-e 2 



|A(t) 



ai\ x 



Api - pi cos(kix + (pi(t) + 9i) sm((kix-(fi(t) + 8i) 

Pi 



(Bi) ai 



2h 



-e 2 



ai\pi sm((kix+<fi(t)+9i 



([E,,B,])« = -A;, 



h 



2e n V 



-A(t) 



Re(A|Q| 2 /2-Qe; 



where pi — \ei\, (pi — arge^ = J dt?p 2 {t') and 6>z = argct/ = const. 
The three second moments are found as 



(ABO* = * 



2e V 



AW |A|q| 2 /2 



-A(t) 



Cov{E u Bi 



-h 



h 



-AW 



Im(A| ei | 2 /2-e ze ;) 



(61) 

(62) 
(63) 

(64) 
(65) 
(66) 



Note also the dumping factors e _A / 2 or e~ A in the expressions of all the above averages. 
The same exponential damping factor e _A occurs in the expressions of mean field energy 
(W{) in time-evolved CS \cxf,t) and time-evolved number states |rij;t). We have 



where 



(Wi) a = ^e" A « [a*(t)a 2 + a^af + Q(t)\ ai \ 2 + 2c,(*)] 
{W l ) n = e(t)e- A ® (l + |) m, 

h 



ai = - 



2e 



. 1 : 



h 

4e 



^(*)N 2 + 



. 1 : 



(67) 
(68) 

(69) 
(70) 



Next we check the Robertson-Schrodinger UR for the one mode fields Ei, Bi and find 
that the obtained three second moments (EU) - fl66l) do minimize it, 



(71) 



h 2 



Thus the time-evolved CS \ot-ut) in nonstationary and/or conductive media are 
Robertson-Schrodinger intelligent states with respect to the photon ladder operator 
quadratures qi, pi, and with respect to the electric and magnetic field components as 
well. In their time development these states can exhibit qi-pi and Ei-Bi correlations and 
squeezing. 
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Concluding Remarks 

In the framework of Choi-Yeon quantization scheme we have constructed ladder operator 
dynamical invariants Aj(t) of electromagnetic field in nonstationary linear media in 
terms of parametric functions that obey the classical nonstationary oscillator 

equation, the modulus \ei(t) \ being subject to the Ermakov equation. The connection of 
ei(t) with the Choi-Yeon [7J[8] parameters 7/, Mi is established. We specified the initial 
conditions ( I49jl under which the eigenstates \azf,t) of Ai(t) (and \ni\t) of A\(t)Ai(t) ) 
represent the time evolution of initial Glauber CS (and photon number states). We 
showed that the medium conductivity and time-dependent permeability act as squeezing 
and correlating factors. For the quantized one mode EMF the first and second moments 
of photon quadrature operators and the electric and magnetic field components in 
\oti\t) are calculated in terms of |e|. It is shown that they are damping in time, and 
the corresponding variances and covariances minimize Robertson-Schrodinger UR. The 
initial photon number states evolve into \ni;t), which are no more eigenstates of the 
(initial) photon number operator. The calculated mean field energy in both \acf,t) and 
\ni;t) exhibits exponential damping factor exp(— A(t)). 

In a similar way one can consider the time evolution in such media of intitial 
canonical squeezed states of electromagnetic field (eigenstates of non-Hermitian 
combination of photon ladder operators) and convince that these states, unlike 
the Glauber CS, remain temporally stable, the varying second moments of photon 
quadrature operators and of the electric and magnetic field components minimizing 
Robertson-Schrodinger UR in all times. 
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